Breached pairing in trapped three-color atomic Fermi gases 
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We introduce an exactly solvable model for trapped three-color atom gases. Applications to a 
cigar-shaped trapped cold fermions reveals a complex structure of breached pairing phases. We find 
two competing superfluid phases at weak and intermediate couplings, each one with two color pair 
condensates, that can be distinguished from density profile measurements. 
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Macroscopic coherent phenomena in matter, such as 
superconductivity and superfluidity, are deep manifesta- 
tions of wave mechanics with consequences not only in 
technology but also for the fundamental understanding 
of vacuum condensates in the standard model. Color 
superconductivity is predicted to occur in quark mat- 
ter at sufficiently high density and low temperatures 
Quarks, having three different colors (red, green, blue) 
and a strong attractive interaction, allow for more di- 
verse pairing patterns compared to the SU(2) Cooper 
pairing in classic metallic superconductors. Such diver- 
sity, likewise, makes it hard to establish the particular 
pairing symmetry favored by nature. With the advent 
of ultracold trapped Fermi g window of opportu- 

nities has opened to address some of these fundamental 
questions, at least, in a qualitative fashion (see for exam- 
ple 01). One can certainly manipulate different atomic 
species and hyperfme states to effectively generate mul- 
ticolor Fermi gases with attractive interactions. 

It is the goal of this paper to investigate the superfluid 
behavior of an imbalanced three-color Fermi gas by means 
of an exactly-solvable pairing model of the Richardson- 
Gaudin (RG) type, derived from the quadratic invariants 
of the SO (6) RG model @, i]. Previous studies using 
standard mean- field density matrix renormalization 
group 0], or Bethe ansatz 0] techniques concentrated 
on the competition between a trionic or barionic phase 
and a color superfluid phase. A main result of our work, 
from the standpoint of collective behavior, is the compe- 
tition between breached-pair (BP) and unbreached-pair 
(UP) superfluid phases in a polarized multicolor Fermi 
gas. As in the two-color case, where density profiles have 
been recently investigated experimentally [8[ and theoret- 
ically , an analogue of the BP or Sarma phase 1(], II | 
appear in multicolor polarized Fermi gases. In this case 
we find a complex structure of breached pairing as well as 
the coexistence of two pair condensates. While the possi- 
bility of coexistence of several superfluid phases has been 
suggested in [l2| using an LDA theory, we predict the ex- 
istence two distinct color fermionic condensates. Within 
our model this is a genuine effect, although care must be 
exercised when contrasted to experiments since interac- 
tions not included in our model could make this phase 



unstable against the formation of a fraction of bound 
trions in the strong coupling limit. However, popula- 
tion imbalance as well as the experimental realization 
of a stable three-color atomic gas with different atomic 
masses and/or different Fcshbach resonances, as recently 
reported in experiments with balanced mixtures of three 
components 6 Li atoms [H, 14 1, could stabilize it. 

Consider the SU(3) color-symmetric Hamiltonian 



H = J2e i N i -gJ2J2 A l A i 



(1) 



for L levels i of energy Si, where cv=(red(i?), green(G), 
and blue(B)) is the color index, TV; = J2 a Nia is the 
number operator of the orbit i, A\ a 

A„„ = (a]J\ , are the pair creators, and g > the pair 
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ing strength. Here atg creates a (canonical) fermionic 
atom in level i with color /3, and e a p 7 is the completely 
antisymmetric tensor in color space. The SO (6) algebra 
(15 generators) is completed by the 9 particle-hole oper- 
ators Ci <a p = a\ a 0np- These 9 operators, which include 
the number operators for the three different colors in the 
level i, N ia = Ci^ aa , close an U(3) subalgebra of S0(6). 

The Hamiltonian |T]) has equal pairing strengths g 
(or scattering lengths), and equal single particle ener- 
gies (masses) for the three colors. It describes a three- 
color Fermi gas with attractive contact interactions in 
the low-density limit. Moreover, the SU(3) symmetry 
is preserved and, thus, the eigenstates are organized in 
degenerate SU(3) multiplets. The SU(3) symmetry, how- 
ever, may be broken by choosing a different combination 
of integrals of motion and, for example, one can generate 
an integrable model of atoms with unequal masses [lj| ■ 

The exact solution of the SO (6) RG model, being an 
algebra of rank three, depends on three sets of spectral 
parameters. The first set includes the usual pair ener- 
gies e a of the S0(6) algebra, while the other two sets, 
composed of the spectral parameters uj a and j a , are as- 
sociated with the SU(3) subalgebra of SO (6). The com- 
plete set of spectral parameters satisfy the generalized 
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Richardson equations 
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The number of spectral parameters in each set is deter- 
mined by the number of particles of each color N a , and by 
the total seniority quantum number v. M — (N — v) /2 
and Q = (Nr +Nr — Ng + v) /2, where we have assumed, 
without loss of generality, > Nr > Nr. The senior- 
ity of level i, Vi, counts the number of unpaired fermions, 
and is defined from Ai a \vi) = 0, JVj \vi) = \fi), 
vi =0,1. The total seniority is v = ~Ylii v i- Hamilto- 
nian (fTJ) preserves the seniority since it can create or de- 
stroy pairs of particles conserving the number parity of 
the level, i.e., for a given configuration each level has an 
even (odd) number of particles where z/, =0 (1). 

The eigenvalues of Hamiltonian |T]) are 



M 



^ = E^ + E 



SiVi 



(3) 



and only depend on the parameters e a . The correspond- 
ing eigenfunctions, though, arc determined by the three 
sets of parameters. The solutions of the Richardson equa- 
tions defines a basis which spans the complete many-body 
Hilbert space of the system. 

Assume a polarized gas composed of TV = Nq + Nr + 
Nr fermionic atoms. In the weak coupling limit the en- 
ergy levels are filled up to the Fermi energy for each color 
SN a , a situation depicted in Fig. 1 A for Nq = 80, 
Nr = 50 and Nr = 20. In this case the seniorities vi are 
equal to 1 for i < 20 and 50 < i < 80, and they are for 
20 < i < 50 and i > 80, defining a clear separation of the 
Hilbert space into regions of odd particle states (i/j = 1) 
and regions of even particle states {yi = 0). When the 
pairing interaction g is switched on, R-B pairs from the 
first region (i < 20) cannot scatter to the second region 
(20 < i < 50) due to Pauli blocking, and they have to 
jump this forbidden region to scatter into the third region 
50 < i < 80. Analogously, G-R pairs of the second region 
have to jump the third forbidden region to scatter into 
the fourth region (i > 80). This configuration, that we 
call BP state, turns out to be the ground state (GS) at 
weak coupling. For larger values of g other configurations 
compete with the BP state. Those configurations, at the 
cost of increasing their kinetic energy, reduce the effect 
of Pauli blocking, therefore, facilitating the pair scatter- 
ing into interior level regions. In particular, the UP state 



depicted in panel B has no blocked interior region, and 
will be the GS of the system at strong coupling. 
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FIG. 1: Occupation number for the BP state (A) and the UP 
state (B). 
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FIG. 2: Phase diagram of a ID partially polarized trapped 
Fermi gas with N = 150, N R = 50, P = {N G - N B )/{N G + 
Nb), and L = 500. 

Although Hamiltonian ((T|) is exactly solvable in any 
dimension, for simplicity we will consider a system of 
N = 150 (Nr = 50) fermionic atoms trapped by a ID 
harmonic potential (of frequency ui) with an energy cutoff 
at E cut = 500huj, implying L = 500 threefold degenerate 
single particle levels. In the weak coupling limit, the 
Richardson equations © decouple into independent sets 
of equations, each one related to the single particle level, 
partially or fully occupied, as discussed above. These 
equations can be solved analytically producing three sets 
of spectral parameters. These parameters are used as an 
initial seed in an iterative procedure in which the cou- 
pling constant g is systematically increased by using the 
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solution of the previous step as the new initial guess. In 
this way the initial solution at weak coupling is evolved 
up to the desire value of g. We performed extensive cal- 
culations to determine the quantum phase diagram of 
this system as a function of the pairing strength g and 
the polarization P = (N G - N B )/(N G + N B ). 

Two color superfluid phases emerge as a function of 
color asymmetry and pairing strength (see fig. [5]) . A first 
order quantum phase transition, due to level crossing, 
separates the BP and UP superfluid phases, which are 
labeled by different sets of seniority quantum numbers. 
On the other hand, there is a smooth crossover between 
the two superfluid phases and a normal, fluctuation- 
dominated state depicted by a thick grey line [l6| . The 
normal Fermi-liquid-like state is dominated by pairing 
fluctuations which arc fully taken into account by the 
exact solution. We adopted the criterium that the nor- 
mal region extends, for a given P, from g = up to the 
value of g for which 10% of the pair energies e a are com- 
plex, meaning that the condensate fraction / is < 0.1 
[171 ] . To study the correlations and structure of the BP 
and the UP states we have chosen two particular points, 
marked with open circles in fig. [5] Notice that within 
our model seniority is a conserved quantum number, dis- 
favoring the formation of bound trionic molecular states 
which arc more likely to appear for very low population 
imbalance and very strong coupling. 
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FIG. 3: Spectral parameters of the BP state with P = 0.8 
and the UP state with P = 0.4, for g = 0.16. 

Figure [3] shows the three sets of spectral parameters 
for these two states. Note that for this value of g = 0.16 
all pair energies are complex with a positive real part, 
implying that the condensate fraction is / = 1, and all 
pairs behave as Cooper resonances as opposed to bound 
molecules of a Bose Einstein condensate (BEC) [l7| ■ In 
both cases the pair energies form two separate arcs in 
the complex plane, indicating the existence of two color 
fermionic pair condensates. In the strong coupling limit 
the two arcs coalesce into a single arc corresponding to 



a single BEC condensate with negative real part of their 
pair energies. The lower arc is overlapping with arcs of 
the two other spectral parameters u a and 7 a which ac- 
count for the couplings in the SU(3) color subspacc. The 
interpretation is that the upper arc of isolated pair en- 
ergies e a describes G-R Cooper pairs, while the lower 
arc corresponds to R-B Cooper pairs. An analysis in 
terms of the eigenvalues of the two-body density matrix 
would lead to a macroscopic eigenvalue of the G-R and 
the R-B pair density matrices and no macroscopic eigen- 
value in the G-B pair density matrix. The appearance of 
these two condensates will be reflected in the occupation 
probabilities that can be calculated using the Hellmann- 
Feynman theorem on the integrals of motion as will be 
explained in a forthcoming paper [iljj ]. 
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FIG. 4: (Color online) Occupation numbers of the BP state 
with P = 0.8 and the UP state with P = 0.4, for g = 0.16. In 
the inset we display the different color contributions to {Ni). 

The occupation probabilities, (Ni), for both states are 
depicted in fig. 01 In the BP state panel we see how the 
G-R pairs avoid the region 50 < i < 90, Pauli blocked by 
the G atoms, to scatter off into the region i > 90. Anal- 
ogously, though less evident due to the smaller number 
of pairs, the 10 R-B pairs avoid the region 10 < i < 50, 
blocked by the R atoms, to scatter off into the region 
i > 50. We may realize that in the latter case the block- 
ing is not perfect, unlike in the two-component gases, be- 
cause of the depletion of the R atoms in the 10 < i < 50 
region due to G-R pairing. The G-B pairing is prevented 
from being realized due to two consecutive blocked re- 
gions. On the contrary, in the UP state pairs do not 
have to avoid blocked regions. The G-R pairs in the 
50 < i < 70 region scatter off into the regions i > 70, 
and R-B pairs in the region i < 30 scatter off into the 
region 30 < i < 50 as well as into the depleted region 
i > 50. These two different physical scenarios manifest 
in the arc geometry of the spectral parameters in fig. [3] 
While the BP state has an inner condensate of 10 R-B 
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and an outer condensate of 40 G-R pairs, the UP state 
has 30 condensed R-B and 20 condensed G-R pairs. 
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FIG. 5: (Color online) Radial density profiles for the BP and 
UP states for g = 0.16. The dash curve corresponds to the 
Thomas-Fermi approximation for the R atoms. The insets 
show the density differences between the 7? and B (dash dark 
gray), and G and R (solid light gray) species. 

An experimental way to uncover the nature of the color 
supcrfluid correlations consists in measuring the density 
clouds of the trapped Fermi gas. Figure [5] shows the ra- 
dial density profiles for these two states normalized to the 
R density as a function of the distance r from the center 
of the trap, in units of the Thomas-Fermi radius R of 
the R species. The insets display the density differences 
between the R and B (dash dark gray), and G and R 
(solid light gray) species. In both cases the outer region 
is dominated by G-R pairing that enforces a maximum 
overlap of the G and R components of the wavefunction. 
However, the BP state shows a richer structure in the 
inner region. In particular, a clear signal of a BP state 
appears as a peak in the density profile of the major- 
ity atoms G, which is magnified in the G and R density 
difference. Two peaks appear in the R and B density 
difference consistent with the decay and revival of the B 
density in regions of R-B pairing. The structure of the 
density differences is smoother in the UP phase due to 
the absence of breached pairing. Regions of constant en- 
ergy differences signal the dominance of the correspond- 
ing pairing phase. In the UP state for 0.0 < r/R < 0.5 
it is dominated by a R-B supcrfluid, while for r/Ro > 1.0 
we have a G-R supcrfluid. 

One can think about measuring the fraction of each 
of the two color pair condensates by exploiting the dif- 
ferences in the color-pair dependent Feshbach resonances 
14j . It would then be possible to use the ramp tech- 
nique as described in [l9(, sweeping the magnetic field 
such that only one class of pairs are transformed into 



bound molecules, allowing for the determination of the 
corresponding fraction of the condensate. A bigger ramp 
would then transform all pairs into molecules, therefore, 
allowing measurement of the complete fraction of the con- 
densate. 
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Si[N iG + N lB + (1 + A)N iR ] - 
^iapCiipa, where A 



